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In a companion paper [ 5uant-ph/9904013| we have investigated several variations of Schwinger's 
proposed mechanism for sonoluminescence. We demonstrated that any realistic version of 
Schwinger's mechanism must depend on extremely rapid (femtosecond) changes in refractive index, 
and discussed ways in which this might be physically plausible. To keep that discussion tractable, the 
^ . technical computations in that paper were limited to the case of a homogeneous dielectric medium. 

' In this paper we investigate the additional complications introduced by finite-volume effects. The 

, basic physical scenario remains the same, but we now deal with finite spherical bubbles, and so must 

' decompose the electromagnetic field into Spherical Harmonics and Bessel functions. We demonstrate 

l/^ ' how to set up the formalism for calculating Bogolubov coefficients in the sudden approximation, and 

, show that we qualitatively retain the results previously obtained using the homogeneous-dielectric 

0^ ' (infinite volume) approximation. 

a^ 

:^ 

§ ■ I. INTRODUCTION 

^3 ■ 

^ '' Sonoluminescence (SL) Q is a phenomenon whose underlying physical mechanism is still highly controversial. In a 
companion paper |Q we have extensively discussed Schwinger's proposed mechanism: a mechanism based on changes 
in the quantum-electrodynamic (QED) vacuum. We discussed Schwinger's original version of the model [||-^, the 
rS Eberlein variant pC|-p^, and proposed a variant of our own [Q. (See also p^-pTj.) 

One of the key features of photon production by a space-dependent and time-dependent refractive index is that 
for a change occurring on a timescale r, the amount of photon production is exponentially suppressed by an amount 
exp(— wr) The importance for SL is that the experimental spectrum is not exponentially suppressed at least out 
to the far ultraviolet. Thus the timescale for change in the refractive index must be of order a femtosecond, and any 
Casimir-based model has to take into account that the change in the refractive index cannot be due just to the change 
in the bubble radius. To achieve this, we adjust basic aspects of the model: We will move away from the original 
Schwinger suggestion, in that it is no longer the collapse from Rmax to Rmin that is important. Instead we postulate 
a rapid (femtosecond) change in refractive index of the gas bubble when it hits the van der Waals hard core 

In this paper we develop the relevant formalism for taking finite volume effects into account. The quantized 
electromagnetic field is decomposed into a set of basis states (spherical harmonics times radial wavefunctions) . The 
radial wavefunctions are piecewise Bessel functions with suitable boundary conditions applied at the surface of the 
bubble. Using these basis states and the sudden approximation we can calculate the Bogolubov coefficients relating 
"in" and "out" vacuum states. In the infinite volume, the discussion of the companion paper is recovered, while 
for physically realistic finite volumes we see significant but not overwhelming modifications. The uncertainties in 
our knowledge of the refractive index as a function of frequency are also significant. We find that we can get 
both qualitative and approximately quantitative agreement with the experimentally observed spectrum, but there is a 
considerable amount of unknown condensed-matter physics associated with the details of the behavior of the refractive 



1 



index. We feel that the theoretical calculations have now been pushed as far as is meaningful given the current state 
of knowledge, and that further progress will depend on looking for the experimental signatures discussed in 

II. BOGOLUBOV COEFFICIENTS 

To estimate the spectrum and efficiency of photon production we decided to study a single pulsation of the bubble. 
We are not concerned with the detailed dynamics of the bubble surface. In analogy with the subtraction procedure 
of the static calculations of Schwinger or of Carlson et al. p^-|20| we shall consider two different configurations. 
An "in" configuration with a bubble of refractive index n™^ in a medium of dielectric constant Coutside, and an 
"out" configuration with a bubble of refractive index n°^l in a medium of dielectric constant toutside- These two 
configurations will correspond to two different bases for the quantization of the field. (For the sake of simplicity we 
take, as Schwinger did, only the electric part of QED, reducing the problem to a scalar electrodynamics). The two 
bases will be related by Bogolubov coefhcients in the usual way. Once we determine these coefficients we easily get the 
number of created particles per mode, and from this the spectrum. (This calculation uses the "sudden approximation" : 
Changes in the refractive index are assumed to be non-adiabatic, see references [p|,p^|-^T| for more discussion.) 

Let us adopt the Schwinger formalism and consider the equations of the electric field in spherical coordinates and 
with a time-independent dielectric constant. (We temporarily set c — 1 for ease of notation, and shall reintroduce 
appropriate factors of the speed of light when needed for clarity.) Then in the asymptotic future and asymptotic past, 
where the refractive index is taken to be time- independent, we are interested in solving 

e{r)d^{d^E)-V^E = Q, (1) 

with e(r) being piecewise constant. We look for solutions of the form 

E = <^{r,t)Yur,m-. (2) 
r 



Then one finds 



e(a2cO) _ (^2$) + hil + 1)$ = 0. (3) 



For both the "in" and "out" solution the field equation in r is given by: 

tdl^ - 9^$ + + 1)$ = 0. (4) 

In both asymptotic regimes (past and future) one has a static situation (a bubble of dielectric ng"^ in the dielectric 
niiquidi or a bubble of dielectric in the dielectric nuquid) so one can in this limit factorize the time and radius 
dependence of the modes: $(r, i) = e*'^*/(r). One gets 

/"+ -!?(/ + 1)^ / = 0. (5) 

This is a well known differential equation. To handle it more easily in a standard way we can cast it as an eigenvalues 
problem 

(^^a + (6) 
where — tuP. With the change of variables / = r^l'^G, so that ^(r, i) ~ e^^*'r^^'^G{r), we get 



1 ^' / 2 



,2 



G +-G + [k^-^]G = 0. (7) 



This is the standard Bessel equation. It admits as solutions the Bessel and Neumann functions of the first type, 
J^{Kr) and N^{Kr), with v = 1 + 1/2. Remember that for those solutions which have to be well-defined at the origin, 
r = 0, regularity implies the absence of the Neumann functions. For both the "in" and the "out" basis we have to 
take into account that the dielectric constant changes at the bubble radius {R). In fact we have 
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^inside " {''^^gasY = dielectric constant of Eir-gas mixture 



(8) 



^outside ~ ''^liquid ~ dielectric constant of ambient liquid (typically water) if r > _R. 



After the change in refractive index, we have 



("gas)^ = dielectric constant of air-gas mixture 



if r < ii, 



(9) 



^outside 



'^'liquid ~ dielectric constant of ambient liquid (typically water) \i r > R. 



In the original version of the Schwinger model it was usual to simplify calculations by using the fact that the 
dielectric constant of air is approximately equal 1 at standard temperature and pressure (STP), and then dealing only 
with the dielectric constant of water {nuquid = \J^outside ~ 1-3). We wish to avoid this temptation on the grounds that 
the sonoluminescent flash is known to occur within 500 picoseconds of the bubble achieving minimum radius. Under 
these conditions the gases trapped in the bubble are close to the absolute maximum density implied by the hard core 
repulsion incorporated into the van der Waals equation of state. Gas densities are approximately one million times 
atmospheric and conditions are nowhere near STP. (For details, see |2l[] page 5437, and the discussion in |^.) For 
this reason we shall explicitly keep track of n™^, and nuquid in the formalism we develop. Defining the "in" and 
"out" frequencies, ujin and ujout respectively, one has 



Gin f in \ 
V \^qasi '^liquid ■, ^in j 



■^in Ain j ( in \ 
"jy ^1/ •Jiyy^qas^inr) 



ifr<R, 



(10) 



^™ [B™ Jv{nuquid Wj„r) + C™ N^iniiqmd Wi„r)] if r > i?. 



Here S™ is an overall normalization. The A™, i?™, and C™ coefficients are determined by the matching conditions 
at R 



'^liquid ^inR), 



(11) 



— Ju'iniiquid ^inR) + N,y' {uuquid ^inR)-, 



(the primes above denote derivatives with respect to r), together with the convention that 

|B|2 + |Cp = l. (12) 

The "out" basis is easily obtained solving the same equations but systematically replacing n™^ by n°"*. There will 
be additional coefficients, S°"*, -B°"*, and C°"*, corresponding to the "out" basis. 

] we adopted a normalization convention such that S°"* was always equal 
to unity. Though that convention is physically equivalent to this one (there are compensating factors in the phase 
space measure), in this paper we find it more convenient to explicitly keep track of this overall prefactor because it 
helps us to compare our finite volume results to the analytically tractable homogeneous model discussed in |^.) 
We now demand the existence of a normalized scalar product such that we can define orthonormal eigenfunctions 

($\$J)=(5*^. (13) 

For the special case Ugas = nnquid, (corresponding to a completely homogeneous space, in which case A = I = B, 
C = 0), it is useful to adopt a slight variant of the ordinary scalar product, defined by 

(01, (/.2) =tn^ I (t)ldo h d^x, (14) 



(In an earlier work jTsI | quant-ph/9805031 



3 



though we shall soon see that this definition will need to be generalized for a time-dependent and position-dependent 
refractive index. (See Appendix A, and compare to the discussion in ^). If we now take the scalar product of two 
eigenf unctions, we expect to obtain a normalization condition which can be written as 

Inserting the explicit form of the $ functions we get 

*k..=«H,„..]) = ^* '^«' ^rnm' u^u, + u,) j vdv J,(nc..r) Unu^T) e»(— .)* (16) 

•J 

— I Ow Oram' n [^i+i^jj \ 6^ ^' (17) 

nu)i nijjj 

= 2 S* Sj 5w Smrn' n S{k., - Kj), (18) 

where we have used the inversion formula for Hankel Integral transforms |^,^ , which can be written as 

' rdr Uk.t) Uk^t) = = 2 ^ ^ ~ = 2 5{n\ n\\ (19) 

this result being valid for Re(y) > — i, and K(i,2) > 0. 

We now compare this to the behaviour of the three-dimensional delta function in momentum space 

53 (k, = ^^''- - ''i^ ^2(^^ _ (20) 

Y^Y:rSK>^,)yira{Q,A,) (21) 
Ira 

5mm', (22) 

to deduce that for homogeneous spaces the most useful normalization is 

2 S* 5^. „ = (23) 

Kill t\j j 

This strongly suggests that even for static but non-homogeneous dielectric configurations it will be advantageous to 
set 

1 (24) 



5{t^^ ■ 


- f^i) 


Hi 


Kj 




- ^3) 




Kj 




- I^j) 


K,i 


Kj 



2n K 

where n is now the refractive index at spatial infinity. 

In our calculations the phase of S is never physically important. If desired it can be fixed by using the well-known 
decomposition of the plane-waves into spherical harmonics and Bessel functions. See, e.g. Jackson |2^ pages 767 and 
740, equations (16.127) and (16.9). The properly normalized plane- wave states of the companion paper are 

[277^/^ V2 w n 

Using the spherical decomposition of plane waves this equals 

E{x, t)^Yl -r exp(-zc.i) ^;™(^(^)' '^(*)) my (26) 

This allows us to identify 

v2 n K 
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thereby fixing the phase, and verifying our normalization from another point of view. 

To confirm that this is still the most appropriate normalization for non-homogeneous dielectrics requires a brief 
digression: The fact that even for homogeneous spaces the usual inner product needs an explicit factor of to get the 
correct normalization is our first signal that the inner product should be somewhat modified for position-dependent 
and time-dependent refractive indices. If we now consider the case of a time-independent but position-dependent 
refractive index, then as explained in Appendix A, the inner product must be generalized to 

(28) 



1,02) = j / £(?■) (/'t do 4>2 d X, 



Taking the scalar product of a pair of eigenfunctions, inserting the explicit form of the $ functions, and recalling that 
K = nuj, we get 



rdr (rii 



rrdr {n,,,f E* A* A, M^rrUA^T'^) 
Jo 



iquid ) >— 



c/ liquid 

^2 "^'^t 



2 {uj, + ujj) E* Ej {B*Bj + C*Cj} {nuqu^d) 

2 + ..,) E* S, {B*B, + dC,}^-^^^^^^ e^(— 
2E*E,{B*B,+C*C,}5{u,-Lu,). 



Jiquid \ 



liquid , liquid 



(29) 

(30) 

(31) 
(32) 



The last few lines follow from some Bessel function identities we have collected into Appendix B. These identities can be 
derived as generalizations of the Hankel integral transform formula (p^, or via use of specific spectral representations 
of the delta function. There are delicate cancellations between surface terms at i?+ and i?^, and the subtle part of 
the calculation involves the surface term at spatial infinity. This calculation is most useful in that it verifies for us 
the normalization condition we need on the "in" and "out" asymptotic states: If we adopt the convention that 



\B\^ + \C\^ = 

then proper normalization of the wavefunctions demands 

1 



1, 



(33) 



(34) 



and we see that it is indeed the refractive index at spatial infinity that is the relevant one for this overall normalization. 

We are now ready to ask what happens if we change the refractive index by making e(r, t) a function of both position 
and time. We are interested in solving the equation 



doie{r,t)doE)^V^E^O, 



(35) 



which is the relevant generalization of (|l|) to the time-dependent case. It should be no surprise that this change affects 
both the conserved density and flux, and that the inner product must again be modified (in what is now a rather 
obvious fashion). See Appendix A for the details. 



2) 



e{r,t) do h d^x, 



St 



The Bogolubov coefficients (relative to this inner product) can now be defined as 



,E^ 



(36) 



(37) 
(38) 
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Where Ej" now denotes an exact solution of the time-dependent equation (|3^) that in the infinite past approaches 
a solution of the static equation (|l]) with e — > ei„(r) and eigen-frequency ujj. Similarly now denotes an exact 
solution of the time-dependent equation (|35| ) that in the infinite future approaches a solution of the static equation 
(|l|) with e —^ eoMt(r) and eigen-frequency uji. The inner product used to define the Bogolubov coefficients has been 
carefully arranged to correspond to a "conserved charge" . With the conventions we have in place the absolute values 
of the Bogolubov coefficient are independent of the choice of time-slice St on which the spatial integral is evaluated. 
With minor modifications, as explained in Appendix A, the inner product can further be generalized to enable it to 
be defined for any arbitrary edgeless achronal spacelike hypersurface, not just the constant time time-slices. (This 
whole formalism is very closely related to the S-matrix formalism of quantum field theories, where the S-matrix relates 
asymptotic "in" and "out" states.) 

Of course, evaluating the Bogolubov coefficients involves solving the exact time-dependent problem (|35|), subject to 
the specified boundary conditions, a task that is in general hopeless. It is at this stage that we shall explicitly invoke 
the sudden approximation by choosing the dielectric constant to be|^ 

e{r, t) = e,„(r) e{-t) + e„„t(r) e{t). (39) 

This is a simple step- function transition from ei„(r) to eout{r) at time t — 0. For t < the exact eigenstates are 
given in terms of the static problem with e = einif), and for t > the exact eigenstates are given in terms of the 
static problem with e = eo„t(r). To evaluate the Bogolubov coefficients in the simplest manner, we chose the spacelike 
hypersurface to be the t = hyperplane. The inner product then reduces to 

(01, (/.2) = « / e(r, i = 0) 0J do 02 d^x, (40) 
Jt=o 

with the relevant eigenmodes being those of the static "in" and "out" problems. (At a fundamental level, this 
formalism is just a slight modification of the standard machinery of the sudden approximation in quantum mechanical 
perturbation theory.) 

There is actually a serious ambiguity hiding here: What value are we to assign to e(r, t — 0)? One particularly simple 
candidate is 

e(r, i = 0) ^ i [ein{r) + eout{r)] = \ [n„(r)2 -f nout(rf \ , (41) 



but this candidate is far from unique. For instance, we could rewrite (39) as 

e(r,i) = cxp (ln{e™(r)} e(-t) + ln{e„„t(r)} e(t)). (42) 
For t 7^ this is identical to (p9[), but for t = this would more naturally lead to the prescription 



e(r,t 0) ^ \J tin(rYout{r) = ?^m(?')?^o«^(^')• (43) 

By making a comparison with the analytic calculation for homogeneous media presented in we shall in fact show 
that this is the correct prescription, but for the meantime will simply adopt the notation 

e(r,i = 0) = 7(n,„(r);no„t(r)) , (44) 

where the only property of 7(^1; n-i) that we really need to use at this stage is that when n\ — ni 

7(n; n) = r? . (45) 

(This property follows automatically from considering the static time-independent case.) 

We are mainly interested in the Bogolubov coefficient /3, since it is |/3p that is linked to the total number of particles 
created. By a direct substitution it is easy to find the expression: 



This was implicit in the earlier work @ [ huant-ph/980503l| . We make it explicit here since this point has caused some 
confusion. 
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(46) 



X / 7 (n,„ (r ) ; Uout [r)) G^"* (n™* , nuqutd , ^out , r) G\? {nll^ , nuquzd , , r) rdr. (47) 

"'0 

(The Srn.-m' arises because of the absence of a relative complex conjugation in the angular integrals for (3. On the 
other hand, the a coefficient will be proportional to Smm'-) To compute the radial integral one needs some ingenuity, 
let us write the equations of motion for two different values of the eigenvalues, ki and K2- 
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If we multiply the first by and the second by G^i we get 

GkiGk2 + ~G^^Gk,2 + — -^{l + ^)^^ GkiGk2 — 0, 

^Ka^Ki + -G^a^^i + ^'^2 ^ + 5)^^ Gi^^G^^ = 0. 

Subtracting the second from the first we then obtain 



Gk2 G^,^Gk.^ 



~ [Gi^^Gk.2 ~ Gi^,^Gki^ + (^2 ~ k1)Gk,^Gi^2 — 0. 



The second term on the left hand side is a pseudo-Wronskian determinant 

W.^.Ar) = Gl^(r)G.,(r) - G',^ir)G.Ar), 



(48) 
(49) 

(50) 
(51) 

(52) 
(53) 



and the first term is its total derivative dW^-^K^/dr. (This is a pseudo-Wronskian, not a true Wronskian, since the 
two functions G^i and G^j correspond to different eigenvalues and so solve different differential equations.) The 
derivatives are all with respect to the variable r. Using this definition we can cast the integral over r of the product 
of two given solutions into a simple form. Generically: 



(k2 - K?) / rdr G^^G^2 



That is 



rdr^^ 
dr 



dr W.; 



dr W^,, 



drW., 



(54) 



(55) 



So the final result is 



rdr Gi^^G^^ = — ^ (^«;if=2 



This expression can be applied (piecewise) in our specific case [equation (^7|)] . We obtain: 



r 7 (n,„ (r) ; Uout (r ) ) G""* (n°a^ , nu.^.d , c^°"* , r) Gj," (n» , , nu,u^d , , r) 

= Trdrj {nil,; n^) GZ^^nlll u;outr)GT {nH^ c^r) 
Jo 

+ / r dr {niiquidf G°J'^{niiquid^^outr)G™{niiquid^^inr) 



(56) 



(57) 



(58) 
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/ V gas I "'gas J 



= R 



7 (n 



gas 1 gas ) 



(59) 



K^^>o«^)2 - (n- 



(u;o„t)2 - (Wj„)2 



(60) 



where we have used the fact that the above forms are weU behaved (and equal to 0) for r = 0. There is an additional 
delta-function contribution, proportional to S^ujin — coout), arising from spatial infinity r = oo. In the case of the P 
Bogolubov coefficient this can quietly be discarded because of the exphcit (uin—tOout) prefactoi]^ For the a Bogolubov 
coefficient we would need to explicitly keep track of this delta-function contribution, since it is ultimately responsible 
for the correct normalization of the eigenmodes if we were to take '^oos- (Here and henceforth we shall 



automatically give the same I value to the "in" and "out" solutions by using the fact that equation (|4^) contains a 
Kronecker delta in I and I' .) Finally the two pseudo-Wronskians above are actually equal (by the junction condition 
(|rT|)). This equality allows to rewrite integral in equation (^7|) in a more compact form 

r drj{n^n{r); Uout [r)) G°"* (n°"* .nuguid, i^out , r) G™ (n^" ^ , niiquid ,ujin,r) 



" " /l*^ A out p 



linTas-^n'^gas) 



1 



l^outV - His '^^n)' i^out)^ - 



— I— 'zn ^out -^Tj ^ 



[{7 Kas;'^gas) " KaD'Knt " {7 K"a- "s-) " i^^asfW 



W[Mn<;,'il LOoutr), JujnTas ^»»r)]j 
[te^Wont)2- (n™ ^,„)2] 



(61) 



(62) 



Inserting this expression into equation ( ^ ) we get 



W[MKal UJoutr), JAKls '^^nr)]l 



[{7 {nils-, '^gas) - ("gas)'}^o«t " {7 "gas) 



Kasrycoi] 



oi(l^out-|-l^in)t 



(63) 



As a consistency check, this expression has the desirable property that /3 as "-gas- That is, if there is no 

change in the refractive index, there is no particle production. We are mainly interested in the square of this coefficient 
summed over / and m. It is in fact this quantity that is linked to the spectrum of the "out" particles present in the 
"in" vacuum, and it is this quantity that is related to the total energy emitted. Including all appropriate dimensional 
factors (c, h) we would have (in a plane wave basis) 



diV(K, 



liquid \ 
jut ) 



73 -t liquid 
" ^out 



liquid ->liquid\\2 i3 -* liquid 



(64) 



Here, since are are interested in the asymptotic behaviour of the photons after they escape from the bubble and move 
to spatial infinity, we have been careful to express the wave-vectors in terms of the refractive index of the ambient 
liquid. This is equivalent to^ 



■^We wish to thank Joshua Feinberg for some insightful questions on this point that caused us to delve into this issue more 
deeply. 

^Remember that when the photons cross the gas-liquid interface their frequency, though not their wave-number, is conserved. 
So we do not need to distinguish ujgas from ujuquid- 
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dN{K, 



liquid\ 
out ) 



liquid 



o/-' liquid -< liquids 
Py^in ^ ^out > 



liquid 



liquid 
out 



liquid 



d J^jri d 



(65) 



If we now convert this to a spherical harmonic basis the angular integrals must be replaced by sums over I, V and 
m,m'. Furthermore we can also replace the dnin and duout by the associated frequencies dujin and dujout to obtain 



dN{uJout) 

dtOmit 



IV mm.' 

In view of our previous definition of the S factors this implies 

dN{uJoiit) 



dujo 



\(3{uJin,UJout)\'^ , 
2~ ^^in^ 



where we have now defined 



\P{uJin,UJout)f = ^ ^ [Plm,l'm'{uJin,(^out)Y 



(66) 



(67) 



(68) 



Im I'm' 



Note that the normalization factors S quietly cancel out of the physically observable number spectrum. Other 
quantities of physical interest are 



N 



duJnu.t 



duJn 



and 



E : 



dN{uJout) 

diOnut 



IjJout duJout- 



Hence we shall concentrate on the computation of: 



(69) 



(70) 



(71) 



Im I'm' 

[{7 {nils', n 



R 



out \ 
gas } 



I out\2\ 2 
\"'gas) S^ou 



{7 {nT„.\n 



out \ 
gas 1 ' "gas ) 



K^s?}^n 



X V(2Z+1) |S,„|2|S„„,|2 



1=1 



H 2 



te^o«t)2-(n» u;,„)2 



(72) 



(Note the symmetry under interchange of "in" and "out" ; moreover / = is excluded since there is no monopole 
radiation for electromagnetism. Also, note that the refractive index of the liquid in which the bubble is embedded 
shows up only indirectly: in the A and S coefficients.) The above is a general result applicable to any dielectric sphere 
that undergoes sudden change in refractive index. However, this expression is far too complex to allow a practical 
analytical resolution of the general case. For the specific case of sonoluminescence, using our variant of the dynamical 
Casimir effect, we shall show that the terms appearing in it can be suitably approximated in such a way as to obtain 
a tractable form that yields useful information about the main predictions of this model. We shall first consider 
the large volume limit, which will allow us to compare this result to Schwinger's calculation, and then develop some 
numerical approximations suitable to estimating the predicted spectra for finite volume. 



III. THE LARGE R LIMIT 



The large R limit is a reasonably good approximation to the physical situation in sonoluminescence, since with 
our new interpretation the radius of the light emitting region is about 500 nm, which is somewhat larger than the 
short distance cutoff on the wavelength {Xmin ~ 200 nm, see ) . Independent of the issue of whether the large R 
limit is a good fit to empirical reality, it is certainly useful in its own right for giving us an analytically tractable 
qualitative understanding of the physics of sudden dielectric changes in large bubbles. (Compare this with the kR ^ 1 
approximation invoked by Eberlein.) 
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If R is very large (but finite in order to avoid infra-red divergences) then the "in" and the "out" modes can both 
be described by ordinary Bessel functions 

G-(n™„ u;, r) = S,„ J,(n™, c.„r/c), (73) 
G°"*(n°:^*, u;, r) - So,* J.(n-'* Uoutr/c). (74) 

We can now compute the Bogolubov coefficients relating these states 

a,,^{Er,E^) (75) 

= S*„ Eout ^'^'"y°"*^ e'("-^— )* <5„, 5„„,, 7 "9-) / ^.«a. io,r.r/c)Mni:l u^outr/c) r dr (76) 

^* I I \ i(u;o„t-t^i„)t r r (in out\ ^{''^^^as^in — ^^^luJout) 
= ^„ ^o„t (Wj„ + Wo„t)e ^ dii' dmm' 7 ("gasi'^gaj " (77) 

= S:„ So„t 7 (n-,; n-'*) + e^("°— <5„, <5(n-, - n™* c.^^*)- (78) 

In terms of a plane-wave basis this is equivalent to 

a(AC„, K^ut) = 7 «a.; '^sa^) f + e'^'^""'"'^-)* _ ^^^^) (79) 

\ "^gas ^gas J 

Note that as — > n™^ the a coefficient has the correct limit: 

aij 2 S*„ Sout (^mm' ^(iJ^in — UJout)- (80) 

Recall that the S factors have been carefully chosen to make sure that the above is simply a three-dimensional delta 
function in momentum space, translated into the spherical-polar basis. That is, in terms of a plane-wave basis the 
large R limit (for rii„ = no„t) is 

a{K.in,K.out) S^{Kin - Kout)- (81) 

The computation for /3 is analogous 

A, = (i?r**,i?r) (82) 

= S,„ 5„„, (^m-^o.O ^.(^„„,+^,„)t ^ / ^»^/c) ^^o^^r/c) r dr (83) 

(84) 



"out (.".'m Li'otit^f^ <Jll' <Jm.-m' I {"'gas : "-gas ) j„ 



= 5o„* 7 <a*) (-^ - -ir) e'(-"+--)* 5™.-™, S{n^:, ^„ - ^o^O- (85) 

\"'gas "'gas/ 

In terms of a plane-wave basis this is equivalent to 

P{n., n,) = 7 ("™.; O (-^ - -ir) e»(-"+-"')* <5^(«,„ + (86) 

\"gas "gas/ 

Comparing this result with the independent calculation of the companion paper [^j finally fixes the overall normal- 
ization of 7(n™,;n°^*): 

7KL;"sa*) =<."sa*- (87) 

Once this normalization is fixed, the large R limit of our finite volume calculation (which we have carried out only in 
the sudden approximation) exactly reproduces the sudden limit of the homogeneous dielectric calculation carried out 
m |]. Readers are referred to that paper for the corresponding numerical estimates of the quantity and spectrum of 
the emitted photons. 

We mention that with considerably more brute force analysis the same results as above can also be obtained directly 
from equation ^73) by formally taking the R oo limit and using the asymptotic formulae for the Bessel functions. 
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To estimate the possible significance of finite volume effects note that in our new version of Schwinger's model we 
have R « Rught-emitung-region ~ 500 nm and take K « 27r/(200 nm) so that KR « Stt 15. (Which is why we 
expect that using the large R limit is a tolerably good approximation, in addition to being much clearer for gaining 
qualitative understanding). To get about one million photons we need, for instance, n™^ « 1 and n°"* « 12, or 
n™, « 2 X lO'' and n™* « 1, or even n™^ « 71 with n™* « 25 @. Note that the estimated values of n°^* are 
extremely sensitive to the precise choice of 7 and the high-frequency cutoff, and that the approximations used in 
taking the large R limit are at this stage uncontrolled. 



IV. BEHAVIOUR FOR FINITE RADIUS: NUMERICAL ANALYSIS 



We now turn to the study of the predictions of the model in the case of finite radius. Unfortunately this cannot 
be done analytically due to the wild behaviour of the pseudo-Wronskian of the Bessel functions. Nevertheless with 
some ingenuity, and a detailed study of the different parts of the Bogolubov coefficient, we are led to some reasonable 
approximations that allow a clear description of the photon spectrum predicted by the model. 



A. The A factor 



The A^, B^, and factors can be obtained by a two step calculation. First one must solve the system (|11|) by 
expressing B and C as functions of A. Then one can fix A by requiring B^ + C'^ = 1 , a condition which comes from 
the asymptotic behaviour of the Bessel functions. Following this procedure, and again suppressing factors of c for 
notational convenience, we find that for the "in" coefficients 



quid 



Bl 



'^liquid 



A 



cr = A 



W[J^{niiquid UJinr), N^{niiquid UJinr)] 
W[J^{niiquid i^inr), J^in'gas UJinr)] 



'^liquid ^in'^ quid 



We are mostly interested in the coefficient A^. This can be simplified by using a well known formula (cf. 
360 formula 9.1.16) for the (true) Wronskian of Bessel functions of the first and second kind. 



(89) 

(90) 
page 

(91) 



In our case, taking into account that for our pseudo-Wronskian the derivatives are with respect to r (not with respect 
to z), one gets for the numerator of A^,: 



W[J^{niiquid U!inr),N^{ni 

iquid '^m^)]i? — '^liquid 



quid 



Hence the A^ factor can be written as 



4/(7r2i?2) 



W[Jv{n'f^as ^inr),Ni,{nUquidOJi7ir)Y + W^[-^i^("-gas ^rnr), Jv(nuquid Wi„r)]2|^ 

Now adopt the notation y = n™^ ujinR/c and yuquid = nuquid i^mR/c = [nuqutd / n''^^^)y . Then 

4/7r2 



Jv{y) N^iyiiquid) 
y Jliy) yuquid Nllyuq^id) 



J Ay) JAyii quid ) 

y J'v{y) yitquid J'^yuquid) 



(92) 



(93) 



(94) 



where in this equation the primes now signify derivatives with respect to the full arguments (y or yuquid)- A similar 
formula holds of course for in terms of x and • Using the standard identities xJ^ (a;) — xJi,-x{x)~vJu{x) and 
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xN[,[x) = xN^-i{x) — vN^[x), applying properties of the determinant, and adopting the notation Min = 'T^/iguid/'T-gasJ 
this can be simphfied to 



4/^2 



Jiy{y) Jvi-N'iny) 

y Jy-l{y) Niny Jv-li-hfiny) 



(95) 



By considering the small argument expansions for the Bessel functions it is relatively easy to see that for small y 
(holding Nin fixed) 



|A™(2/^0,M„)|2^(M„)2- + O(y). 



(96) 



On the other hand, for large values of the argument y the asymptotic forms of the Bessel functions can be used to 
demonstrate that 



iAr(y->oo,M„)|2 



'^^gas "^liquid 



(97) 



Numerical plots of [A^p show that it is an oscillating function of y which rapidly reaches this asymptotic form. The 
mean value for large arguments is simply: 



1 /■^'^ 



^"'gas "'liquid 



= 1. 



(98) 



While these results are general, for the particular application to SL that we have in mind it is the small y behaviour 
that is most relevant. Also, keep in mind that this large y asymptotic formula holds for y large but ignoring dispersive 
effects (that is, assuming a frequency independent index of refraction). If we model dispersive effects by a Schwinger- 
like cutoff where the refractive index drops to unity (see below) then above the cutoff we will have = 1 holding as 
an identity. 



B. The Pseudo— Wronskian 



Use the simplified notation in which x — uJoutR/ y = fig^s ^inR/ c. In these dimensionless quantities, after 
making explicit the dependence on R and c, and inserting the particular choice of 7 motivated by the large- i? limit, 
equation (^2|) takes the form: 



/ in 2 I out 2 
.2 1^ |2 1^ |2 / "'gas *^ ~ '^gas U 



^out ^in 



I '—'in I I ^out I 



nTas a; + y 



F{x,y). 



(99) 



Here F{x, y) is shorthand for the fimction 



F{x,y)=J2i^l + l) \An' \Ai 



Ju{x) Ju(y) 
x Jl{x) y Jliy) 



1=1 



(100) 



where in this equation the primes now signify derivatives with respect to the full arguments (x or y) . It is convenient 
to define a dimensionless Bogolubov coefhcient, and a dimensionless spectrum, by taking 



i?2 



mx,y)\^ 



(101) 



so that 



dN{x) 
dx 



A j^in rv-fOut 
^ "'gas"' gas JQ 



dy 



^in I ^ I '^out I ^ 



(102) 



The total number of photons is then 
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N 



A :j^in ^ out 
^ "•gas"'gas JO 



dx 



dy 



The total energy emitted is given by a very similar formulc 

he 1 



E 



^^"'gas ^ "'gas'^gas JO 



dx 



dy X -ri: 



(103) 



(104) 



In order to proceed in our analysis we need now to perform the summation over angular momentum. Although the 
infinite sum is analytically intractable, we can easily demonstrate that it is convergent and can physically argue that 
the lowest angular momentum modes will dominate the sum. Consider the large order expansion [v ^ x dX fixed x) 
of the Bessel functions. In this limit one gets psl : 



Jy{x) 



This can be used to obtain the asymptotic form of the pseudo-Wronskian appearing in equation (100) 



W,{x,y) 



Ju{x) Ju{y) 

X Jlix) y Jl{y) 

Jv{x) J„{y) 
X J^+i{x) y J„+i{y) 

{x^ - y'^) ( xy 



27r(z^)l/2(j,+ 1)3/2 



(f) 



2v+l 



(105) 

(106) 
(107) 
(108) 



where we have used the standard recursion relation for the Bessel functions Jl,{z) = i'J^{z) — zJ^^i(z). This indicates 
that the sum over v is convergent: the terms for which [xy/v'^) < 1 are suppressed. Whatever the values of x and y 
are, for sufficiently large angular momenta this asymptotic form guarantees the convergence of the sum over angular 
momenta. 

Everything so far has been predicated on the absence of dispersion: the refractive index is independent of frequency. 
In real physical materials the refractive index is known to fall to unity at high enough energies. (Sufficiently high 
energy photons "see" a vacuum inhabited by effectively-free isolated charged particles. The manner in which the 
refractive index approaches unity is governed by the plasma frequency, and the location of this physical cutoff is 
governed by the resonances present in the atomic structure of the atoms.) This situation is far too complex to be 
modelled in detail, but it is easy to see that an upper bound on emitted photon energies implies an upper bound on 
the allowed angular momentum modes: Basically, if one supposes the photons to be produced inside or at most on 
the surface of the light emitting region, then the upper limit for the angular momentum (as seen at spatial infinity) 
will be attained by photons emitted tangentially from the edge of the light emitting region: this maximal angular 
momentum is the product of the radius of the light emitting region times the maximum observed "out" momentum. 
Then one gets: 



outside 



'hserved , 



X R 



= RK, 



observed' 



(109) 



For sonoluminescence Kobserved is of order 27r/(200 nm). Since the light emitting region is known to be approximately 
500 nm wide we shall be most interested in the case KR « 57r « 15, with a corresponding maximum angular 
momentum Imax approximately 15. Under these conditions, the bulk of the radiation will be into the lowest allowed 
angular momentum modes. The precise value of the angular momentum cutoff Imax is sensitive to the details of both 
the frequency cutoff in refractive index, and the size of the light emitting region. For instance, in some of Schwinger's 
papers he took K k, 27r/(400 nm) in which case (taking again R « 400 nm) Imax ~ 57r/2 « 7. Whatever ones views 
as to the precise value of this cutoff it is clear that the emitted radiation is limited to low angular momenta. 

A subtlety is that this is the angular momentum as measured at spatial infinity (in the ambient liquid — water). This 
is not the same as the angular momentum the photons have while they are inside the bubble (since it is frequency, not 
wavenumber, that is conserved when photons cross a timelike interface [spacelike normal] ).p| Taking this into account 



*For a flash occurring at minimum radius he/ R ~ 0.4 eV. 

^Contrast this to a spacelike interface (timelike normal; sudden temporal change in the refractive index) for which it is the 
wavenumber, not the frequency, that is conserved across the interface. During photon production we are dealing with a spacelike 
interface, whereas when the photons escape from the gas bubble we are dealing with a timelike interface. 
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jinside gQ^ joutside jj ^ -i c /'1in\ 

''max ~ ''max ~ tXJXobserved ~ iO- V-"^-"^"; 

"^liquid ^liquid ^liquid 

We now discuss how to take this observationaUy based cutoff in angular momentum and translate it into a statement 
about the cutoff in the refractive index. 



C. Implementation of the cutoff 

If we adopt a Schwinger-like momentum-space cutoff in the refractive index, then because we have defined the 
variables x and y partly in terms of the refractive index, we must carefully assess the meaning of these variables. In 
terms of momenta, Schwinger's cutoff is 

riinin) = riin Q{Kin - n) + I Q{k - Kin), (HI) 

nout{n) = nout Q{Kout - k) + 1 6(k - Kout), (112) 
This implies that the photon dispersion relation uj{k) has a kink at k = K , and that we can write 

io^n{^^) = Q{K^n " «) + + c(k - i^.„) Q{n - K,n), (113) 

\ ^171 / 

uJout{K) = Q{Kout - k) + + c{k - K,^t) Q{k - Kout). (114) 

'^out \ ^out J 

Finally, the variables x and y generalize (actually, simplify) to 

x = K.outR/c; y = K.inRlc, (115) 

so that 

niniy) = riin Q{y* - y) + 1 Q{y - ?/*), (116) 

n-outix) = Uout Qix* - .t) + 1 6(a; - a;,), (117) 

where = KoutR/c; y* = KinR/c. Now all these changes do not affect F(x,y), which is why we defined it the 
way we did, but they do affect the prefactors appearing in equation (p9|). An immediate consequence is that the {x, y) 
plane naturally separates into four regions and that |/3(x,?/)p — in the region x > x^ and y > y*. We shall soon 
see that the two "tail" regions (x < x^,y > y,) and (x > x^,y < y^) are relatively uninteresting, and that the bulk 
of the contribution to the emission spectrum comes from the region (x < x^,y < y^,)^ 

In the infinite volume limit this is an exact statement, since in that limit one can show (see |^ and the discussion 
below) that F{x,y) — *■ G{x) S{x — y) so that the support of the spectral integral is exactly the line segment x = y 
with X — y < min{a;,, ?/*}. 

Finally, when it comes to choosing specific values for a;* and y*, we use the fact that the variables x and y are 
related to the angular momentum cutoff discussed in the previous subsection to set 

= = 15. (118) 

IT'liquid 



If one is too enthusiastic about adopting the sudden approximation then the integral over these tail regions will be divergent. 
This, however, is not a physical divergence, but is instead a purely mathematical artifact of taking the sudden approximation 
all the way out to infinite frequency. The integral over these two tail regions is in fact cut off by the fact that for high enough 
frequency the sudden approximation breaks down. As a practical matter we have found that the numerical contribution from 
these tail regions are small. 
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D. Working along the diagonal 



To study in more detail the behaviour of the function F{x, y) when higher angular momentum modes are retained 
one can perform a Taylor expansion of F{x^ y) around x = y. 



lim 



— lim 



Ju{x) J^{y) 
X Jl{x) y Jl{y) 



y {x — y) x^y [x — y) 

(x) + {x- y)Jl{x) 
X Jlix) X Jlix) + {x- y)[Jl{x) + X J'^ix)] 



lim 

x^y 



{x - y) 

Mx) -Kix) 

X Jl{x) Jl{x)+x J'J{x) 
= J,{x)[JUx) + X J'J{x)] - X Jlixf. 

The derivatives can be eliminated by using the well known recursion relations. 



W^ix,y) 

lim — — = Ji,(x) 

x^y (X — y) 



'{y^-x^y 




— X 




X 





-Ji,{x) - J^+i{x) 



2vJ^{x)J^+i{x) - X [Jl{x) + Jl^-^ix)] 



(119) 

(120) 

(121) 
(122) 

(123) 
(124) 



For sake of simplicity we shall use an equivalent form of equation (124) where lower order Bessel functions appear 

WAx.y) 



hm 

x^y [X — y) 



2vJ^{x)J^^i{x)~x [Jl{x) + Jl_^{x)\ 



(125) 



This result shows that, as expected, each term of F{x, x) is finite along the diagonal and equal to zero at x = y = 0. 
Moreover 



D{x) = F{x,x) =^(2/ + l) 



1=1 



\{2l + l)Ji+i/2{x)Ji_i/2{x)-X Jf+i,2{x)+ Jf_^/2{x) } 



(126) 



This sum can easily be checked to be convergent for fixed x. [Use equation (105).] With a little more work it can be 
shown that 



lim D{x) = -i^. 



The truncated function obtained after summation over the first few terms (say the first ten or so terms) is a long and 
messy combination of trigonometric functions that can however be easily plotted and approximated in the range of 
interest. A semi- analytical study led us to the approximate form of D{x) 



D{x) 



1 



27r2 250 



(127) 



A confrontation between the two curves in the range of interest is given in the figure below. 
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FIG. 1. Plot of the exact D{x) against its approximated form (dotted curve) in the range < x < 14 



E. The factorization approximation 



To numerically perform the integrals needed to do obtain the spectrum it is useful to note the approximate factor- 
ization property 



F{x,y) 



X + y x + y 



G 



y 



(128) 



That is: to a good approximation F{x, y) is given by its value along the nearest part of the diagonal, multiplied by 
a universal function of the distance away from the diagonal. A little experimental curve fitting is actually enough to 
show that to a good approximation 



F{x,y)^D 



V 



sm' 



^(3[x-y]/4) 



(3[x-y]/4)2 



(129) 



From the plot we show below it is easy to check that the function F{x, y) is quite well approximated by this factorized 
form. We feel important to stress that this is approximation is based on numerical experimentation, and is not an 
analytically-driven approximation. (In the infinite volume case we know that F{x,y) (constant) x S{x — y). The 
effect of finite volume is effectively to "smear out" the delta function. In this regard, it is interesting to observe that 
the combination sin^(x)/(7rx^) is one of the standard approximations to the delta function.) Our approximation is 
quite good everywhere except for values of x and y near the origin (less than 1) where the contribution of the function 
to the integral is very small. 
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FIG. 2. Transverse fit: An orthogonal slice of F{x,y) intersecting the diagonal at {x,y) = (3,3). Here F(3 + z,3 — z) is 
plotted in comparison with [sin^(3a/2)]/(3a/2)^. The solid line is the function, and the dotted line the analytic approximation 



F. The spectrum: numerical evaluation 

We have now transformed the function F{x,y) into an easy to handle product of two functions 

^ 27r2 16000 + (a; + y)6 {3[x-y]/4)^ ' ^ ' 

We exhibit tridimensional graphs for both the exact (apart from the approximation of truncating the sum at a 
finite /) and approximate forms of the function F{x,y). We have chosen the case of R = 500 nm (corresponding to 

= 15 n°g'^l/nuquid as previously explained). 
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FIG. 3. Plot of the exact F{x,y) in the range < a; < 12, < j/ < 12. The jagged behaviour along the diagonal is a 
numerical artifact, as the function is known to be smooth there. 




FIG. 4. Plot of the approximated F{x, y) in the range < x < 12, < j/ < 12 



The dimensionless spectrum, based on equations (|72| ) and (|99|), is 

2 



dx 



„in 2 I out 2 
"oas ~r "nas f 



^ ^gas^gas JO \ ^gas ^ + ^gas V 



x + y\ sin^(3[x^y]/4) 



(131) 



where n°"*(a;) and n™^g{y) are now the appropriate functions of x and y (See equations (117) and ( |116| )). We have 
also manually inserted a factor 2 to account for the photon polarizations. 
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As a consistency check, the infinite volume hmit is equivalent to making the formal replacements 



sin2(3[a; - y]/4) Air 



(3[.T - y]/4)2 3 

and 



S{x-y), (132) 



The first replacement can be formally justified as follows. It is known that a sequence of smooth functions approxi- 
mating the delta function is given by 



indeed, one get 



/.(X) = ^ (134) 



lim f,ix)^dix). (135) 



Then, it is straightforward to show that 

sin^{3[x~y]/A) ^ sin^jR ^KasU^^n - ng^*^out]/(4c)) 
{3[x - 2/]/4)2 {R 3[nf,,cu„, - n-*u;o„t]/(4c))2 

^ ^ K^Kas^^n - <,"*C.„„i]/(4c)) 

= fs{x-y). (136) 

Doing so, equation (131) reduces to the spectrum obtained for homogeneous dielectrics in j^. Indeed 

dN 1 (n™ ~n-*)2 



dx 37r n™,n°Jf* 



x^ e(x, - x). (137) 



With these consistency checks out of the way, it is now possible to perform the integral with respect to y to estimate 
the spectrum for finite volume, and similarly to perform appropriate double integrals with respect to x and y to 
estimate both total photon production and average photon energy. In our companion paper we showed that in the 
infinite volume limit there were two continuous branches of values for n™^ and n""* that led to approximately one 
million emitted photons with an average photon energy of 3/4 the cutoff energy. If we now place the same values of 



refractive index into the formula ( |131| ) derived above, numerical integration again yields approximately one million 
photons with an average photon energy of 3/4 times the cutoff energy. The total number of photons is changed by at 
worst a few percent, while the average photon energy is almost unaffected. (Some specific sample values are reported 
in Table I.) The basic result is this: as expected Q, finite volume effects do not greatly modify the results estimated 
by using the infinite volume limit. Note that hn^nax is approximately 4 eV, so that average photon energy in this 
crude model is about 3 eV. 



in 
"'(jas 


OJlf 

"•(jas 


Number of photons 


{E)/fin,nax 


2 X 10* 


1 


1.06 X 10*^ 


0.803 


71 


25 


1.00 X 10'^ 


0.750 


68 


34 


1.06 X 10" 


0.751 


9 


25 


0.955 X 10" 


0.750 


1 


12 


0.98 X 10" 


0.765 



Table I: Some typical cases. 
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In addition, for the specific case 71^"^ = 2 x 10^, n™* = 1, we have calculated and plotted the form of the spectrum. 
We find that the major result of including finite volume effects is to smear out the otherwise sharp cutoff coming from 
Schwinger's step-function model for the refractive index. Other choices of refractive index lead to qualitatively similar 
spectra. These results are in reasonable agreement (given the simplicity of the present model) with experimental data. 




100000 
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12 



14 



FIG. 5. Spectrum dN/dx obtained by integrating the approximated Bogolubov coefficient. We integrate from y = Q to 
y, — 11.5 and plot tlie resulting spectrum from a; = to x = 14.5. For Uout = 1 and R = 500nm the relation between the 
non-dimensional quantitiy x and the frequency v is x u ■ 10.5 ■ 10~^^s. So x ^ 11.5 corresponds to v ^ 1.1 PHz. The curve 
with the sharp cutoff is the infinite volume approximation. Finite volume effects tend to smear out the sharp discontinuity, but 
do not greatly affect the total number of photons emitted. 



V. DISCUSSION 



The present paper presents calculations of the Bogolubov coefficients relating the two QED vacuum states appro- 
priate to changes in the refractive index of a dielectric bubble. We have verified by explicit computation that photons 
are produced by rapid changes in the refractive index, and are in agreement with Schwinger in that QED vacuum 
effects remain a viable candidate for explaining SL. However, some details of the particular model considered in the 
present paper are somewhat different from that originally envisaged by Schwinger. Based largely on the fact that 
efficient photon production requires timescales of the order of femtoseconds we were led to consider rapid changes in 
the refractive index as the gas bubble bounces off the van der Waals hard core. It is important to realize that the 
speed of sound in the gas bubble can become relativistic at this stage ^ . 

We feel that theoretical computations along these lines have now been pushed as far as is meaningful given the cur- 
rent experimental situation. We have now constrained Casimir-like mechanisms for sonoluminescence into a relatively 
small region of parameter space, and it is experiment, rather than theory, that is likely to lead to further advances. 
We argue, both here and elsewhere, that Casimir-like mechanisms are viable, that they make both qualitative p^ ] 
and quantitative predictions, and that they are now sufficiently well defined to be experimentally falsifiable. (Possible 
extensions of the model will require much more detailed condensed matter information, such as experimental data 
regarding the actual refractive index inside the bubble as a function of time, space, and frequency.) 
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In conclusion, the present calculation (limited though it may be) represents an important advance: There now can 
be no doubt that changes in the refractive index of the gas inside the bubble lead to production of real photons — the 
controversial issues now move to quantitative ones of precise fitting of the observed experimental data. We are hopeful 
that more detailed models and data fitting will provide better explanations of the details of the SL effect. 
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APPENDIX A: GENERALIZING THE INNER PRODUCT 

For the differential equation 

e doidoE) ~ V^E = 0, (Al) 
it is a standard exercise to write down a density and flux, 

p = e{EldtE2-E2dtE*), (A2) 



j^ElVE2-E2VEl, (A3) 
and to then show that, by virtue of the differential equation ( [Al[ ), these quantities satisfy a continuity equation 

dtp-V-j = 0. (A4) 



Suppose now one has two solutions of the differential equation (Al), one can then define an inner product 

{Ei,E2) = ~i eJ^{E*dtE2 ~ E2dtEl) , (A5) 

where the integral is taken over a constant-time spacelike hypersurface. By virtue of the above, this inner product is 
independent of the time t at which it is evaluated. 

Now what happens if the dielectric is allowed to depend on both space and time? First the differential equation of 
interest is generalized to 

do{e{x, t)doE) ~ V^E = 0. (A6) 

Second, the density and flux become, 

p = e(r, t) {EldtE2 - E2dtEl) , (A7) 



J = EIVE2 - E2\7El . (A8) 

By virtue of the differential equation (|A6|), 

dtp = El dt{eix,t)dtE2) - E2 dMx,t)dtEl) (A9) 

= EIV^E2 - E2V^El (AlO) 

= V • {EIVE2 - E2S/EI) (All) 

= V-j. (A12) 
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Which imphes that these generalized quantities satisfy a continuity equation 

dtp-V-j = 0. (A13) 



This imphes that the generahzed inner-product [for two solutions Ei and E2 of the equation (A6) for time-dependent 
and space-dependent dielectric constants] must be 

{Ei,E2)^-'-J^e{x,t){EldtE2-E2dtEl). (A14) 

By the continuity equation this inner product is independent of the time t at which the integral is evaluated [provided 
of course, that Ei and E2 both satisfy (A6)]. This construction can be made completely relativistic. Define a 
four-vector Jf^ by 

J^^ip;f). (A15) 
Then for any edgeless achronal spacelike hypersurface E there is a conserved inner product 

{Ei,E2) = -i J e{x,t) J'' dJ:^. (A16) 



APPENDIX B: SOME BESSEL FUNCTION IDENTITIES 

The inversion formula for Hankel Integral transforms is [ |2^ , p3[ | 

' r dr Mmr) Mn2r) = = 2 1^^!!^ = 2 5(^1 4), (Bl) 

^Kl K2 [Kl + K2) 

this result being valid for Re{y) > — i, and A(i^2) > 0- On the other hand, another well-known standard result is 

Kin Ki-i 



(See, for instance (24.88) on page 142 of Spiegel |2^. Cf. also |24 , formula 11.3.29 on page 484. This is of course a 
special case of the more general pseudo-Wronskian analysis presented in the text.) This is enough to imply 

R [kiJI{kiR)Ji,{k.2R) - K2Jy{niR)Jl(K2R)] _ &{ki - K2) 
R^oo K2 — k\ K2 

By using the asymptotic forms of the Bessel functions (see, for example, equations (24.103)-(24.104) of Spiegel [p6[) 
this is equivalent to the two spectral identities 

sin(fci?) 

hm — ^ — -^■n5{k), (B4) 

and 

hm 2^ ^ 0. (B5) 

R-^00 k 

These spectral identities, together with the known asymptotic forms of the Bessel functions, then let us generalize 
(B3) above to obtain both 

R[KiNl{KiR)N,iK2R) - K2N,{kiR)NI{k2R)] S{ki-K2) 

lim 5 5 = ■ (Bd) 

-R^oo K2 — nf y/Ki K2 



and 



R[,liJliKiR)N,{n2R) ~ H2JA^ilR)K{>i2R)] „ 

hm ^ 5 = (J. (B7j 



These are the key equations needed to complete the argument leading to the correct normalization of the static 
eigenmodes. [See equation (^.] 



22 



t 
t 

t 

§ 

[1] 

[2] 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 



[14] 

[15] 
[16] 
[17] 

[18] 
[19] 
[20] 
[21] 
[22] 

[23] 
[24] 
[25] 
[26] 



E-mail: liberati@sissa.it 
E-mail: visser@kiwi.wustl.edu 
E-mail: belgiorno@mi.infn.it 
E-mail: sciama@sissa.it 

B.P. Barber, R.A. Hiller, R. Lofstedt, S.J. Putterman Phys. Rep. 281, 65-143 (1997). 

S. Liberati, M. Visser, F. Belgiorno, and D.W. Sciama, Sonoluminescence as a QED vacuum effect. I: Physical Scenario, 



quant-ph/99040i; 



Schwinger, Proc. Nat. Acad 
Schwinger, Proc. Nat. Acad 
Schwinger, Proc. Nat. Acad 
Schwinger, Proc. Nat. Acad 
Schwinger, Proc. Nat. Acad. Sci 
Schwinger, Proc. Nat. Acad. Sci 
Schwinger, Proc. Nat. Acad. Sci 



Sci, 
Sci, 
Sci, 
Sci, 



89, 4091-4093 (1992). 

89, 11118-11120 (1992). 

90, 958-959 (1993). 
90, 2105-2106 (1993). 
90, 4505-4507 (1993). 

90, 7285-7287 (1993). 

91, 6473-6475 (1994). 



Eberlein, Sonoluminescence as quantum vacuum radiation, Phys. Rev. Lett. 76, 3842 (1996). quant- ph 9506023 

Eberlein, Theory of quantum radiation observed as sonoluminescence, Phys. Rev. A 53, 2772 (1996). quant-ph 9506024 
Eberlein, Sonoluminescence as quantum vacuum radiation (replies to comments), Phys. Rev. Lett. 77, 4961 (1996); 78, 
2269 (1997). |quant-ph/9610034| 
C.S. Unnikrishnan and S. Mukhopadhyay, Phys. Rev. Lett. 77, 4960 (1996). 
A. Lambrecht, M-T. Jaekel, and S. Raynaud, Phys. Rev. Lett. 78, 2267 (1997). 
N. Garcia and A.P. Levanyuk, Phys. Rev. Lett. 78, 2268 (1997). 

S. Liberati, M. Visser, F. Belgiorno, an d D.W. Sciama, Sonoluminescence: Bogolubov coefficients for the QED vacuum of 



a collapsing bubble, quant-ph/9805023, 

S. Liberati, M. Visser, F. Belgiorno, and D.W. Sciama, Sonoluminescence as a QED vacuum effe ct, quant-ph/9805031 . 
S. Liberati, M. Visser, F. Belgiorno, and D.W. Sciama, Sonoluminescence and the QED vacuum, quant-ph/9904008 
F. Belgiorno, S. Lib erati, M. Visser, and D.W. Sciama, Sonoluminescence: Two-photon correlations as a test for thermality. 



quant-ph/99040iq . 

C. E. Carlson, C. Molina-Pan's, J. Perez-Mercader, and M. Visser, Phys. Lett. B 395, 76-82 (199 7). hep-th/9609195 
C. E. Carlson, C. Molina-Pan's, J. Perez-Mercader, and M. Vi sser, Phys. Rev D56, 1262 (1997). [hep-th/9702007|. 



C. Molina-Pan's and M. Visser, Phys. Rev. D56, 6629 (1997). |hep-th/9707073 
J. Holzfuss, M. Riiggeberg, and A. Billo, Phys. Rev. Lett. 81, 5434 (1998). 

H. Bateman, Higher Transcendental Functions, Vol II, (McGraw-Hill, New York, 1953); See especially equation (8.1.1) 
page 5. 

J.D. Jackson, Classical Electrodynamics, (Wiley, New York, 1975); see equation (3.112) on page 110. 

M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, (Dover, New York, 1965). 

A. Jeffrey, Handbook of Mathematical Formulas and Integrals, page 219, (Academic Press, San Diego, 1995). 

M. Spiegel, Mathematical handbook of formulas and tables, (McGraw-Hill, New York, 1968). 



23 



